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Chapter 1 

Electron transport in one dimension 



This chapter gives a simple introduction to scattering problems in one dimension and their 
relation to the conductance of one-dimensional systems at the level of introductory quantum 
mechanics. Purely one-dimensional systems are of limited practical use, but their simplicity 
allows one to introduce the basic physics using simple mathematics only. Sec. 1.1 recapitulates 
the concept of a quantum current, as defined in introductory quantum mechanics. The Landauer 
formula expresses the electronic conductance or resistance of a system as a scattering problem. 1 
A naive derivation is given in Sec. 1.2. 2 Sec. 1.3 introduces some elementary folklore of 
scattering theory, the transfer matrix and the scattering matrix, using the ubiquitous square 
barrier potential as an example. Techniques to solve scattering problems for any one-dimensional 
potential are introduced in Sec. 1.4. These will be generalized to three-dimensional problems in 
the next chapter. In particular I focus on so-called "mode matching" , which is the most accessible 
technique. I discuss its relation to the "Green function" and the "embedding" techniques, 
two other currently used methods. The techniques are illustrated by a simple model, which 
consists of a straightforward discretization of the one-dimensional Hamiltonian. 3 This model is 
mathematically equivalent to a simple tight-binding model of a quantum wire of atoms, as is 
shown in Sec. 1.5. 4 

1.1 Quantum currents 

1.1.1 Probability currents 
The probability current J(x, t) 



x It is also called the Landauer- Biittiker formalism. The first idea to relate the conductance of a small system 
to its scattering matrix comes from Landauer. Biittiker made a generalization to multi-terminal devices, which 
have more than two electrodes. 

2 Secs. 1.1 and 1.2 are adaptations from my lecture notes for the introductory bachelor course in quantum 
mechanics. 

3 The model is didactical and not designed to be the most efficient way to solve this scattering problem. The 
math is easy, however, and the general ideas can be applied to more complicated systems. 
4 Secs. 1.3-1.5 consist of elements of lecture notes I am writing for a master course. 
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is defined such that 



dt 



J(a,t) - J(b, t) 
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describes the change in the probability P a b(t) of finding a particle in the region a < x < b at 
time t, where 



Pab(t)= / \*(x,t)\ Z dx, (1.3) 
J a 

provided the wave function ^(x,t) describing the particle is normalized. Proving this is trivial: 
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Going from the first to the second line one uses the Schrodinger equation ih "^ x '^ 



h 2 d 2 V{x t) 
dt 2m &P h 

V(x)ty(x,t) and its complex conjugate (the potential terms cancel), and from the second to the 
third line one uses partial integration. Setting b = a + dx with dx infinitesimal, allows one to 
write Eq. 1.2 as 



dp(x,t) dJ(x,t) 



dt 



dx 



(1.4) 



with p(x,t) = \^f(x,t)\ the probability density. You might recognize Eq. 1.4 as a continuity 
equation, which describes the relation between a density and a current. 

Probability currents may seem rather abstract, but they are easily related to something 
more familiar. Suppose the particle has a charge q, then the expected charge found in the region 
a < x < b at time t is Q a b{t) = qP a b(t)- 5 Defining the electrical current as I(x, t) = qJ(x, t), 
Eq. 1.2 can be rewritten as 



dQabjt) 

dt 



I(a,t)-I(b,t). 



(1.5) 



This makes sense; the rate of change of charge is given by the difference between the current 
flowing in from one side minus the current flowing out from the other side. 

A nice thing is that, even if the wave function cannot be normalized, like the wave function of 
a free particle, the probability current according to Eq. 1.1 is still a well-defined quantity. Free 
particles often enter in scattering problems, where we are interested in quantities like reflection 
and transmission coefficients. Since the latter can be directly defined in terms of probability 
currents, we can get away with using non-normalizable wave functions. 6 



1.1.2 Stationary states 

Suppose now that ^(x,t) describes a stationary state, i.e. 



&(x,t) = ip(x)e 



(1.6) 



5 This is an expectation value in the quantum mechanical sense. One starts the wave at time and at time t 
one measures whether the particle is in the region a < x < b. By repeating this "experiment" over and over, one 
can calculate the probability P a b(t)- Q a b(t) is then the average charge found in this region from these repeated 
"experiments" . If you have problems imagining this then think of a particle emitter that sends out a pulse of 
many (independent) particles. The averaging is done automatically and the average charge is what you measure. 

6 For quantum purists: one can work with normalizable wave packets. The math then becomes ugly, and in 
the proper limit the physical results will be the same as when using plane waves. 
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Then one finds from Eq. 1.3 

dPabjt) _ 

dt 

and from Eqs. 1.2 and 1.1 

J(x,t) = J = constant. (1.7) 

The probability current is constant, i.e. independent of position and time. 
For example, consider a free particle with the wave function 

ip(x) = A e ikx . (1.8) 

From Eq. 1.3 we calculate 

P ab = \A\ 2 (b-a). (1.9) 

Since P a b is the probability of finding the particle in the interval between x = a and x = b, i.e. 
an interval of length b — a, we can interpret \A\ 2 as the probability density per unit length. It 
is also called the particle density 7 

p=\A\ 2 . (1.10) 
The probability current is easily calculated from its definition, Eq. 1.1 

J= hk 2 = hk p 

m m 



(1.12) 



According to de Broglie's relation p = hk is the momentum of the particle and 

hk p 
m m 

is then the velocity of the particle. The electrical current is given by 

I = qJ = qvp, (1-13) 

which is the usual definition of an electrical current, namely charge x velocity x density. For the 
wave function of Eq. 1.8 both velocity and density are constant, so the wave function describes 
a uniform current. Suppose q > 0; then if k > the current flows to the right, if k < 0, the 
current flows to the left. From now on we assume that k > 0. 
Now let's go to the more complicated wave function 

V»(x) = Ae ikx + Be~ ikx , (1.14) 

with A, B constants. The associated probability current is 

J = — \A\ 2 - — \B\ 2 , (1.15) 
m m 

which is interpreted as a right going current minus a left going current. In a scattering problem 
one would interpret the first term on the right handside of Eq. 1.14 as the incident wave and 
the second term as the reflected wave. Eq. 1.15 is then interpreted as the difference between 
incident and reflected currents 

J = Jin~ Jr. (1.16) 
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Figure 1.1: Cartoon representing a general one-dimenssional scattering problem. In the left 
region the potential is a constant V(x) = Vl, in the middle region the potential V(x) can be 
anything, and in the right region the potential is a constant V(x) = Vr. The middle region is 
called the scattering region. The left and right regions are called the left and right leads. 
In the left lead we have an incoming wave Ae LX and a reflected wave Be~ lkLX and in the right 
lead we have a transmitted wave Fe lkRX . 



The reflection coefficient R is defined as the ratio between reflected and incident currents 

R=^ = K- (1-17) 

Jin \A\ 2 

Consider the scattering problem shown in Fig. 1.1. In the left region we assume that the 
potential is constant V(x) = Vl, in the middle region the potential can have any shape V(x), 
and in the right region the potential is again constant V{x) = Vr. The solution in the left region 
is given by Eq. 1.14 with k replaced by k^ 

The solution in the right region is given by the transmitted wave 

ij)(x) = Fe lkRX ; x in right region, (1-19) 

with 



kR= V^W) (120) 

One can calculate the transmitted current as 

j T = ^\ F \\ (1.21) 
m 

The transmission coefficient T is defined as the ratio between transmitted and incident currents 

2 



Jt = vr\FT 

Jin V L \Af 



T =^ L = z^rT2> ( L22 ) 



7 Using a beam of iV particles in which each particle is independent of the others and is described by the wave 
function of Eq. 1.8, the particle density is N \A\ 2 , which is the number of particles to be found per unit length. 
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using Eq. 1.12. 

From the fact that the current has to be independent of position everywhere, see Eq. 1.7, it 
follows from Eqs. 1.16 and 1.21 that 



This relation expresses the conservation of current, or: "current in = current out" (reflected 
plus transmitted). No matter how weird the potential in the middle region is, the current going 
into it has to be equal to the total current coming out of it. No particles magically appear or 
disappear in the middle region. From the definitions of Eqs. 1.17 and 1.22 it is shown that Eq. 
1.23 is equivalent to 



i.e. the reflection and transmission coefficients add up to 1. Since these coefficients denote the 
probabilities that a particle is reflected or transmitted, this simply states that particles are either 
reflected or transmitted. 

1.2 Quantum conductance 
1.2.1 Tunnel junction 

The device shown in Fig. 1.2 is called a tunnel junction. The left and right regions consist of 
metals and the middle region consists of an insulator material, usually a metal-oxide. 8 Such 
devices can be made in a very controlled way with the middle region having a thickness of a 
few nm. One is interested in electrical currents, i.e. the transport of electrons through such 
junctions, or more generally in the current-voltage characteristics of such a device. 9 On this 
small, nanometer length scale electrons have to be considered as waves and quantum tunneling 
is important. Nano-electronics is the general name of the field where one designs and studies 
special devices that make use of this electron wave behavior. 

We start with the simplest possible one-dimensional model of a tunnel junction. The atoms 
of a material attract electrons by their nuclear Coulomb potential. The electrons in low lying 
energy levels are localized around the atomic nuclei and form the atomic cores. The atomic 
valence electrons experience a much weaker effective atomic core potential, which is the sum 
of attractive nuclear and repulsive core electron terms. If the atoms are closely packed and 
the material is sufficiently simple, all these atomic potentials add up to a total potential that 
is relatively constant in space. As a qualitative level the potential for electrons in a material 
can be approximated by a constant, which is what we are going to do in the following. 10 The 
constant potential depends on the sort of atoms a material is composed of, so it is different for 
every material. The potential in the tunnel junction of Fig. 1.2 along the transport direction 
can then be represented by a square barrier, as shown in Fig. 1.3. 

8 Scanning tunneling microscopy (STM) uses a tunnel junction between the probe tip and a surface, where the 
middle region is simply vacuum. 

9 The device is called MIM, which stands for metal-insulator-metal. Using magnetic metals the device can be 
applied as a magnetic field sensor, or in MRAMs (magnetic random access memories). 

10 This approximation is often used for simple metals such as the alkali's, aluminium, silver and gold. The 
constant potential approximation is also called the jellium approximation. It does not hold for complicated 
metals such as the transition metals or for covalently bonded materials, such as silicon or carbon. To be fair, it 
doesn't even hold very well for simple metals if one is interested in quantitative results. 




Jr + Jt- 



(1.23) 



1 = R + T, 



(1.24) 
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Figure 1.2: Schematic representation of a tunnel junction. The yellow balls represent atoms 
of a metal, the blue balls represent atoms of an insulator. The left and right regions stretch 
macroscopically far into the left and right, respectively. The electron waves in the metal are 
reflected or transmitted by the insulator in the middle region 
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Figure 1.3: Simple approximation of the potential along the transport direction of a tunnel 
junction, see Fig. 1.2. In the metal (left and right regions) the potential is constant, V(x) = V\. 
In the insulator the potential is also constant, V(x) = Vq, where Vq > V±. The incoming, 
reflected and transmitted waves are given by Ae ikx ,Be- ikx and Fe ikx . 
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Figure 1.4: The potential when a bias voltage U is applied between the left and right leads. This 
changes the potential of the right region by AV = —eU to V\ — AV with respect to the left 
region; see Fig. 1.3. The voltage drop is indicated schematically. If the bias voltage is small, i.e. 
AV <C Vo — Vi, then we can still use the transmission coefficient T calculated for the unbiased 
square barrier (given by the dashed line). 



1.2.2 The Landauer formula 

According to Eq. 1.22 the transmitted electrical current is given by 

It = hnT, (1.25) 

using the definition I = qj (the charge q of an electron is — e). The incoming current Ii n is 
given by Eq. 1.13 as 

Iin = -evp. (1-26) 

How large are the velocity v and density p of the incoming electrons? To answer that we must 
ask the more basic question: how is the incoming current in a device created? In an experimental 
setup this is done by applying a voltage difference U between the left and right regions. The left 
and right regions are metals, which can be connected to the two ends of a battery, for instance. 
This results in a potential drop AV = —ell between left and right regions, as shown in Fig. 
1.4. We suppose that the temperature is zero and the metals are non-magnetic, so we have spin 
degeneracy. Then vp is given by the simple expression 

Tin 7rn 

This expression is derived in the next section. 

If the potential drop AV in Fig. 1.4 is small compared to the barrier height Vo — Vi, we can 
use the unbiased square barrier potential from Fig. 1.3 to calculated the transmission coefficient 
T. 11 This is the so-called linear response regime. Eqs. 1.25-1.27 then give for the transmitted 
electrical current, also called the tunneling current 

Ir = f h UT, (1.28) 



11 One can calculate the transmission coefficient T for the tilted barrier of Fig. 1.4, but this calculation is more 
complicated, and for small the result is almost the same as for the square barrier. 
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which is a remarkably simple expression! If we define the conductance Q as current divided 
by voltage, we get 

It e 2 

g = w = ^ < L29 > 

Since T is just a dimensionless number between and 1, e 2 /7rh has the dimension of conductance. 
It is the fundamental quantum of conductance; its value is e 2 /irh 7.75 x 10~ 5 Q . 

Eq. 1.29 is called the Landauer formula; it plays a central role in nano-electronics. 13 As 
it stands here, it is valid for a one-dimensional, spin degenerate system at low voltage and at 
not too high a temperature, but it can be generalized. 

1.2.3 Simple derivation of the Landauer formula 



I will give a very simple derivation of Eq. 1.27. We have to do a little bit of solid state physics, 
but I use only simple introductory quantum mechanics language. Spin degeneracy means that 
each energy level can be filled with two electrons. The non spin degenerate case is relevant for 
magnetic materials. I let you work out that case yourself. 



The Pauli exclusion principle and the Fermi energy 

The left and right regions of a tunnel junction consist of metal wires, see Fig. 1.2. These wires 
are supposed to be very, very long compared to the size of the middle region. In a simple-minded 
model the potential of a metal wire looks like Fig. 1.5. The potential is approximately constant 
inside the wire and it has steps at the beginning and end of the wire to keep the electrons in. 
The energy levels of this square well potential are, 14 

*2l,2 2 _2 fi 2 

The spacing between the energy levels, E n — E n -\, scales as 1/L 2 with the length L of the 
wire. If L is large, the spacing becomes very small, so from a distance the energy level spectrum 
almost looks like a continuum, as illustrated by Fig. 1.5. 
The wave functions are given by 



sin k n x = — )= (e iknX - e - iknX ^j . (1.31) 



These are not exactly what we need, because they correspond to standing waves, whereas we 
need traveling waves to describe currents, see Eq. 1.8. For the incoming current we only need 
the exp(ifc n x) part. Setting A = l/(iy/2L), the corresponding electron density according to Eq. 
1.10 is 

p=\A\ 2 = ± L . (1.32) 



If you are more used to working with resistances, the resistance 1Z is the inverse of the conductance, i.e. 
1Z — 1/Q, so the quantum of resistance is 7vh/e 2 ~ 12.9kfi. 
13 R. Landauer, Philosophical Magazine 21, 863 (1970). 

14 This is actually the solution for an infinitely deep square well, whereas you might think that we need the 
solution for a finite square well. However if the well is very wide and not too shallow, the infinite well is an 
extremely good approximation. 
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Figure 1.5: Schematic drawing of the potential and the energy levels of a long wire. The points 
— L/2 and L/2 mark the beginning and the end of the wire. The spacing between the energy 
levels is so small that the energy spectrum almost looks like a continuum. Ep marks the Fermi 
energy, i.e. the highest level that is occupied in the ground state by an electron. 



The wire is full of electrons since each of the atoms in the wire brings at least one electron 
with it. Filling the energy levels according to the Pauli principle, and having N electrons in 
total, the highest occupied level is En. The highest occupied level in the ground state is called 

2 

the Fermi energy or Ep. 



Incoming and tunnel currents 

Go back to the tunnel junction and fill its potential profile, Fig. 1.3, with electrons from the 
left and right wires. This is shown in Fig. 1.6. The Fermi energy Ep in the left and right 
regions is the same. The exclusion principle then tells us that there can be no flow of current. 
Any electron on the left side that would try to go to the right side finds an energy level that is 
already occupied by an electron, which excludes any other electron from going there, and vice 
versa. This confirms what we know from everyday life; in an unbiased system no current flows. 16 
Now apply a bias voltage between left and right regions as in Fig. 1.4. The result is shown 
in Fig. 1.7. The bias voltage lowers the right region in energy. Suddenly the electrons on the 
left side that occupy energy levels E n with Ep — AV < E n < Ep find empty levels with that 
energy on the right side. They can tunnel through the barrier to occupy these levels. Provided 
the potential drop AV is small, we can approximate the transmission coefficients of all these 
electrons by T at an energy E = E F . 17 The incoming current of Eq. 1.26 has contributions 
from aU electrons with energies between Ep — AV and Ep. 

I in = -2ep y n- ( L33 ) 

E F -AV<E n <E F 

The factor of 2 is there because there are two electrons in each level. 



In an ordinary metal the Fermi energy is of order 10 eV with respect to the lowest energy level of the valence 
electrons. 

16 One can also reverse the argument. If the Fermi energies on the left and right side would be different, then 
a current would flow. However this current would be short-lived. By sending electrons from left to right one 
occupies a level on the right side, and de-occupies a level on the left side. This would go on until the highest 
occupied levels on the left and right are the same; in other words until an equilibrium is reached. 

17 Again, this is valid in the linear response regime. 
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where 



Figure 1.6: Tunnel junction where left and right regions are filled with electrons. The Fermi 
energies Ep on the left and right side are identical. The exclusion principle forbids electrons to 
trespass from left to right or vice versa. 

This sum in Eq. 1.33 is rather awkward, but by a trick we can turning it into an integral 

EVn = -V^ n y = - VVA/c 
7T ^— ' L 7T z — ' 

« - J vdk, (1.34) 

Ak = k n - k n -i = ^, (1.35) 

see Eq. 1.30. Turning the sum into an integral is allowed because L is very large, so Ak is 
tiny. The lower and upper bound of the integral in Eq. 1.34 should correspond to the energies 
Ep — AV and Ep, whereas the integral is over dk, which is again awkward. We can however 
turn it into an integral over dE, using the following trick 

v = — = — - ^) - 1 dE (136) 
m h dk h dk' 

In Eq. 1.34 it gives 

/ vdk = -^-dk = / Ue 

JE F -AV J Ep—AV " "-k JE F -AV " 

= ^{Ep-{E F -AV)) = ^AV. (1.37) 

Collecting Eqs. 1.37, 1.34 and Eq. 1.32 in Eq. 1.33, we find for the incoming current 

Ay 



(1.38) 

This is the required expression for the incoming current, see Eqs. 1.26 and 1.27. The tunnel 
current is then given by 

eAV 

I T = I in T = —T, (1.39) 

TTtl 

where the transmission coefficient T needs to be calculated for the energy E = Ep. Note that 
with AV = —ell, where U is the potential difference (in Volts), this corresponds to Eq. 1.28. 
The Landauer formula, Eq. 1.29, is then derived straightforwardly. 
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Figure 1.7: Tunnel junction with an applied bias voltage. All the levels occupied by electrons in 
the left region with an energy Ep — AV < E n < Ep correspond to empty energy levels in the 
right region. The electrons in these levels can tunnel from though the barrier from left to right. 

1.3 Transfer and scattering matrices 

The Landauer formula expresses the conductance in terms of a transmission coefficient. In other 
words, the problem of finding the conductance becomes a problem of solving the scattering 
problem. In this section I review some basic elements of scattering theory in one dimension. 
You probably know what a scattering matrix is, but read through the section anyway, if only to 
refresh your memory. 

1.3.1 Transmission of a rectangular barrier 

Again I explain the concepts using elementary quantum mechanics only and the rectangular 
barrier as an example. Consider Fig. 1.8 and let the middle region run from x = — a to x = a. 
We are trying to find a solution to the time-independent Schrodinger equation 



0. 



2m dx 2 

Since the potential is piecewise constant, we can write the solution as 



Ae ikx + Be- ikx ; 
Cef + De-i*; 
Fe ikx + Ge- ikx ; 



x < — a 
-a < x < a 
x > a 



(1.40) 



(1.41) 



with 



k 



2mE 



2m(V Q - E) 



n 



(1.42) 



For < E < Vq, both k and r\ are real numbers. 18 Provided we choose the constants A-G such, 
that the function ip(x) is continuous and differentiable everywhere, it will be a solution to the 
Schrodinger equation for all x, including the boundaries x = ±o. Continuity of ip(x) at x = — a 
gives the relation 



Ae -ika + Be ika 



Ce- va + Dei", 



18 Most of the equations actually also hold for < Vo < E (scattering over the barrier) or for Vq < < E 
(scattering from a potential well). In those cases r\ will be a purely imaginary number. 
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Figure 1.8: Scattering from a rextangular barrier. The transfer matrix M relates the coefficients 
of the waves at the right to the waves at the left of the barrier. 

whereas continuity of ^"0( x ) at ar = —a gives 

ikAe- ika - ikBe ika = ^Ce^ - r)De^ a : 
These two equations can be combined into the matrix equation 

4 ' MW ( C n ) (1.43) 



B J \ D 



ik+ri \ (ik-n)a ( ik-rj \ (ik+y)a 

2ik J e ^ 2ik ) e 

ik-r) \ (ik+n)a ( ik+y \ -(ik-rj)a 

2ik I e l 2ik j e 



In a similar way, continuity of tp(x) and ^ip{x) at x = a gives the matrix equation 

(M mw(f) (l44) 



D J \ G f 

where the matrix M^ 2 ) can be obtained from M^ 1 ) by replacing —a by a, ik by r] and rj by ik. 
Combining Eqs. 1.43 and 1.44 gives an equation of the form 

1 ^ M ( F n \ (1.45) 



B ) \ G t 

M = M (1) M (2) , 

where the 2x2 matrix M is the product of the two matrices of Eqs. 1.43 and 1.44. 

In a scattering problem we set our boundary conditions by hand, or, if you wish, by choosing 
specific experimental conditions. We define an incoming wave Ae lkx , with A a fixed value that 
gives the density of particles corresponding to our experimental source, see Eq. 1.10. At the 
same time we set G = 0, which means that we have no incoming wave Ge~ lkx from the right 
(experimentally this is easy, even for a theoretician). Comparing Eqs. 1.22 and 1.45 we find for 
the transmission coefficient 



T= — 

vl 



1 



(1.46) 
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Figure 1.9: Approximating a barrier of any shape by a series of steps. The transfer matrix M 
is given by M = MiMi • • • M N . 

Since the potential is symmetric, we actually have vr = V£ here, but we let the more general 
form of Eq. 1.46 stand for later. It is a simple exercise to work out the matrix element M\\ and 
show that one obtains the usual textbook expression for the transmission coefficient 



r 



,2 



!-l--[^ L )sinh 2 (2^) 



-i 



(1.47) 



Comparing Eqs. 1.17 and 1.45 gives the reflection coefficient as 



R=P^. (1.48) 
1.3.2 Transmission of a general barrier 

The algorithm formulated in the previous section can be extended in order to describe the 
transmission through a barrier of a more general shape. The matrix M as in Eq. 1.45 is called 
the transfer matrix. It relates the modes in the right region to the modes in the left region. 
Once we know the transfer matrix, we can calculate the reflection and transmission coefficients, 
see Eqs. 1.48 and 1.46. For a single step in the potential, the transfer matrices are given by 
Eqs. 1.43 and 1.44 for a step up and a step down, respectively. Since the square barrier can be 
viewed as a step up (at x = —a), followed by step down (at x = a), its transfer matrix is simply 
a product of the transfer matrices of the individual steps, as is expressed by Eq. 1.45. The idea 
can be applied to a potential of any shape, approximating the potential by a series of steps, as 
is illustrated in Fig. 1.9. 

The transfer matrix M is given by the product of the transfer matrices of each step 

M = MiMi • • • Mjv, (1.49) 

where each Mj can be calculated by defining 



* = ^ 2 "'"^'>- E ', (1.50) 

and using a properly modified Eq. 1.43 or Eq. 1.44. I leave the details up to you. In essence this 
is the transfer matrix algorithm. It is used quite frequently in numerical calculations solving 
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scattering problems in various branches of physics; quantum mechanics, optics, acoustics, etc.. 
Obviously it is best suited for layered materials, i.e. systems that consist of layers of different 
materials stacked on top of one another. If the r/j get large, which is for high barriers and/or 
low energies, one might get into numerical problems because of the exponentials in Eqs. 1.43 
and 1.44. There are some special tricks to handle these, but I won't go into details here. For 
systems described in terms of real atoms, the techniques discussed in the following sections are 
better suited. 



1.3.3 The scattering or S matrix 



The transfer matrix M is closely related to an elementary concept in scattering theory, called 
the scattering matrix. Whereas the transfer matrix relates the modes on the right to the modes 
on the left, the scattering matrix relates outgoing modes to incoming modes. For the example 
shown in Fig. 1.8, the modes that are coming into the barrier are Ae lkx (from the left) and 
Q e -ikx (f rom the right) and the modes that are going out from the barrier are Be~ lkx (to the 
left) and Fe lkx (to the right). The scattering matrix is defined by the relation 



B 
F 



S' 



A 
G 



(1.51) 



Comparison to Eq. 1.45 gives for its matrix elements 

M 2i 



ii 



Mis 



Mn' ~ 22 Mn 

I have put a "prime" on S' because the scattering matrix S as it ordinarily used is defined as 

fur 



(1.52) 



o of o I of 

oil = on, 012 = 4/ — O12, 

V V R 

c / c' . c _ c' 

021 — 4/ — o 21 , O22 — o 22 . 

V V L 



(1.53) 



The reason for introducing the ^Jv factors is because one wants the scattering matrix S to 
reflect a basic conservation law, namely the conservation of current. In a stationary problem 
the current going out must be the same as the current coming in, since otherwise one would get 
an accumulation or depletion of particles, as discussed Sec. 1.1. 



J out 
2 



v L \B\* + v R \F 



Jin ^ 

VL |^ 2 + ^|G| 2 



( V^ B \ 


2 


( V^A \ 






\ V^G J 



From Eqs. 1.51 and 1.53 it is easily shown that this holds only if 

sts = ss f = 1. 



(1.54) 



(1.55) 



In other words, conservation of current means that the scattering matrix is unitary. 

It is custom to write the scattering matrix as 



r t' 
t r' 



(1.56) 
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potential barrier 




Figure 1.10: Discretization of a potential barrier. An equidistant grid is applied with A = 
Xi + \ — Xi. The potential is different from a constant for x\ < x < xm- The coupling to the left 
and right regions (the leads) can be included by extending the grid by two points, which defines 
the boundary zones, indicated in gray. 

From Eqs. 1.46, 1.48, 1.52 and 1.53 it follows that the transmission and reflection coefficients 
are given by 



R 



(1.57) 



hence the name transmission and reflection amplitudes for t and r, respectively. Note that the 
relation 1 = R + T, see Eq. 1.24, then simply reflects the unitarity of the scattering matrix. 



1.4 Mode matching 

Although the transfer matrix algorithm is quite general, it becomes quite tedious if the potential 
in the leads is not constant, which is when we want to make an atomistic model of a tunnel 
junction, see Fig. 1.2. In this section I will discuss a technique that is easier to generalize, called 
mode matching. By means of an introduction, I will discuss this technique in its simplest 
form and take Fig. 1.10 as starting point. The potential is discretized as before, but the grid is 
extended by one point into the left and right regions. 

The basic idea is to discretize the whole Schrodinger equation, Eq. 1.40, including the kinetic 
energy. Approximating the second derivative by a simple, first order finite difference one obtains 



2m 



A 2 



Vrfi = 0, 



(1.58) 



where and Vi are shorthand notations for ip{xi) and V(xi), and A = xi + \ — Xj. I only consider 
equidistant grids here. As is usual in a scattering problem, i runs from — oo to oo, so we have 
an infinite number of equations. This is not only awkward, but also unnecessary. The potential 
is localized in space, i.e. V(x) differs from a constant only for x\ < x < X]y. For the left and 
right regions, x < x\ and x > xn, we already know the solutions to the Schrodinger equation. 
They are simple plane waves, as indicated in Fig. 1.10, with 



^2m{E - V L ) ^/2m(E - V R ) 
k L = ; k R = . 



1.59) 
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We need a way of matching these "modes" to the wave function in the region of the potential 
barrier (the scattering region). 

Let's start with xo, and put the origin there, so xo = 0. The finite difference Schrodinger 
equation for i = is 

E^ + ^2 {^i - 2 ^o + 4>-i}- Vo^o = 0- (1-60) 

We know that for x < the wave function has the form ip(x) = Ae lkljX + Be~ lkLX , so 

V>_! = Ae~ ikLA + Be ikLA = Ae~ ikLA + (ip - A) e ikLA , (1.61) 

where the last result follows from the fact that the wave function has to be continuous at x = 0. 
Remember that in a scattering problem we assume that we know the incoming wave, so A is 
fixed. Eq. 1.60 can then be rewritten as 

+ 2^{^~ ^° + ~ = 2^ A { &lkLA ~ C ~ lkLA } ■ (L62) 

The term on the left handside now only contains ip i with i > 0, and the terms at the right 
handside can be considered as the source of the incoming wave. This takes care of the left 
boundary. 

Now focus upon the right boundary; for i = N + 1 we have 

H 2 

E ^N+l + 7^^2 {^ N + 2 ~ 2 ^n+i + V'at} ~ V N+1 tp N+1 = 0. (1.63) 



We now use 



to get 



E1>n + i + ^2 {^N + ie tkRA - ^N+i + ^n] - V N+ i^ N+1 = 0. (1.65) 

Note that we have again used our knowledge of the scattering boundary conditions. We have 
assumed only a transmitted wave in Eq. 1.64 (no incoming wave from this side) and wave 
function continuity at the boundary. The term on the left handside of Eq. 1.65 only contains 
ipi with i < N + 1. This takes care of the right boundary. 

For i = 1,... ,N we have no problems and we can use Eq. 1.58. We can collect these 
equations together with Eqs. 1.62 and 1.65 and summarize the problem as 

(£I-H')V> = q. (1.66) 

Here ip is a vector containing the coefficients ipf, i = 0, . . . , N + 1, q is the "source" vector of 
length N + 1, where all the coefficients are zero, except the first one 

and H' is a (N + 2) x (N + 2) Hamiltonian matrix. All off-diagonal matrix elements are zero, 
except the sub- and super-diagonal ones 
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All diagonal matrix elements are identical to that of the original finite difference Hamiltonian 

%,i = -^& + V *> (1-69) 
except the first and the last one, which are modified to 

H' Qfi = -J^ + V + X L (E), 
H' N+1 , N+1 = -^ + V N+1 + E R (E). (1.70) 

with 

The quantities T, L / R (E) are called the self-energies of the left and right leads. 19 They take 
care of the proper coupling of the potential barrier to the outer regions, and contain all the 
information we require about the outer regions. The self-energy depends upon the energy of the 
incoming and scattered waves, see Eq. 1.59. Note that, while Eq. 1.58 represents an infinite 
dimensional problem, by introducing the self-energy (and the source), we have reduced it to a 
finite, N + 2 dimensional problem, Eq. 1.66. That can be solved using standard algorithms 
for solving linear equations. 20 

Once you have solved Eq. 1.66, the only thing remaining is to extract the transmission and 
reflection amplitudes. The transmission amplitude is simply given by the wave function at the 
right side of the barrier, normalized to the incoming wave, and normalized with the velocities 
(to attain a unitary scattering matrix, see Eq. 1.53) 

(1.72) 

The reflection amplitude is similarly determined from the wave function on the left side minus 
the incoming wave, normalized to the incoming wave 

r=^Z*. d.73) 

Some care should be taken in determining the velocities. Since we have discretized the 
Schrodinger equation, it is consistent to discretize the expression for the current, Eq. 1.1, in a 
similar way 




^(^^-^k-^)- (L74) 

This should actually give a position independent result, since we are considering a stationary 
problem. For a simple plane wave Ae lkx this expression gives 



_ ih\A\ 2 f ikA ikA \ 
J ~2mA\ e 6 J' 



19 a somewhat unfortunate term descending from Green function folklore. 
20 In this case it is particularly simple, since the matrix is tridiagonal. 
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and comparison to Eqs. 1.11 and 1.12 gives 

ih 

V ~ 2mA 



Note that the continuum limit liniA^o of this expression gives Eq. 1.12. 
The source term, Eq. 1.67, can then be simplified to 



Qo = -£~ V L- (1-76) 



In addition, from Eqs. 1.75 and 1.71 one can relate the velocity to the self-energy 

2t 



2A 

v L/R = -—ImZ L/R (E). (1.77) 



1.4.1 Green function expressions 

Using Green functions the mode matching results can be put into a very compact, albeit some- 
what obscure, form. Define a matrix 

G(E)=(EI-H'y 1 . (1.78) 

It is called a Green-function matrix. Note that it has a dimension N + 2, as has the modified 
Hamiltonian matrix H'. One can also define the infinite dimensional retarded Green- function 
matrix related to the original infinite dimensional Hamiltonian 

G r (E)=[(E + i7 ] )l-H}- 1 , (1.79) 

where r\ is (real, positive) infinitesimal. The advanced Green-function matrix is defined as 21 

G a (E) = [G r (E)]K (1.80) 

For z a complex number in the lower half plane, the matrix elements of G and G r in the 
scattering region are identical. 

G i j(z) = G r iJ (z);i,j = 0,...,N + l (1.81) 

A proof of this can be found in the literature. Note that the modified Hamiltonian matrix H' 
is non-Hermitian, essentially because the self-energy S is not real, see Eqs. 1.70 and 1.71. One 
can show that the eigenvalues of H' are not real and lie in the upper half complex plane. Thus, 
G(E) is a well-defined quantity for real energies E (unlike G r ). It has the retarded boundary 
condition build into it and one does not need the +ir/ trick. 22 
The definition of G allows us to write, 

V'jv+i = G N+lfi (E)q , (1.82) 



21 Since E is an eigenvalue of H, EI — H is singular and its inverse does not exist for real energies E. 
Adding/subtracting an imaginary it] avoids the singularities. In textbooks on scattering theory it is shown that 
this has a physical meaning. G r can be used to construct the retarded wave function, which consists of a wave 
coming in from the source to the target, and waves scattered out from the target. This is the physical solution. 
G a gives the reverse, i.e. waves scattered into the target, plus waves going into the source. This is unphysical, 
but can sometimes be useful in formal mathematical manipulations. After having constructed the wave function, 
the linit lim^o can be taken, so r/ only serves as an intermediate to reflect the boundary conditions. 

22 Mathematicians would call G(E) an analytical continuation of G r (E). All the poles of G(E) are in the upper 
half plane, which makes it analytical on the real axis, and a retarded form. 
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see Eq. 1.66. Eqs. 1.72, 1.76 and 1.82 then lead to a compact expression for the transmission 
amplitude 

t= l -^Jv#uZG N+1 , (E). (1.83) 

An expression of this type is called a Fisher-Lee expression. It relates matrix elements of the 
scattering matrix to matrix elements of the Green function. 

For Green function junkies, we can even make it fancier and write, using Eq. 1.77 

t = 2iv/-ImS« G N+lfi (E) V~ImS L . (1.84) 

This allows the transmission coefficient to be written as 



T = t*t = 4(ImS ii )G^ +li0 (ImS L )Gg )iV+1 (1.85) 

with all quantities evaluated at the fixed energy E. This expression is known as the Caroli 
expression. 

Introducing a Green function to tackle this simple one-dimensional problem is like using a 
sledgehammer to crack a nut. Green functions expressions can however be modified to include 
multiple dimensions, a large bias voltage, interactions with vibrations and/or between electrons, 
etc., where they give compact (though not necessarily practical) expressions. In case of a large 
bias voltage, the Caroli expression is also known as the NEGF expression, where NEGF stands 
for non-equilibrium Green function. 23 



1.4.2 The embedding potential 

If you prefer differential equations over linear algebra, we can take the continuum limit of Eq. 
1.66, i.e. limA^o- It is tricky, but straightforward to obtain 

mix) + £ ( %1 + 8{x) d Jm _ S (x - a) d -^P-) (1.86) 
2m \ ax z ax ax J 

-V(x)ip(x) - {6(x)X L (E) + 6(x - a)Z R {E)} ij>(x) = ihv L A6(x), 

where I have put the left boundary at x = and the right boundary at x = a. The S's have the 
form 

X L/R (E) = (1.87) 

The term on the right handside of Eq. 1.86 describes the source term, which only exists at the 
left boundary. On the left handside there are ^-function terms that only exists at the boundaries 
of the barrier with the outside regions. They take care of the coupling to the outside region. 
The ^ terms in the kinetic energy ensure a continuous derivative across the boundaries. The 
potential term between { } is called the embedding potential. The formalism is known in the 
literature as the embedding formalism. 24 The differential equation needs to be solved for ip(x), 
but only over the finite domain < x < a. The coupling to the outside regions is taken care 
of by the boundary terms. Eqs. 1.72 and 1.73 give the transmission and reflection amplitudes 
with ip{a) replacing V'tv+i an d ip{0) replacing ip Q . 



23 This is in contrast to the linear response regime, where the Green functions are ordinary, i.e. equilibrium 
Green functions. 

24 Be careful however, as in different fields the phrase "embedding" is attached to different things. 



20 



CHAPTER 1. ELECTRON TRANSPORT IN ONE DIMENSION 



Alternatively, one may use Green function expressions also in this case. Define the Green 
function G(x, x' , E) as the solution of the equation 

(E-H'^G(x,x',E) = 6(x-x'), (1.88) 

where H' is the operator of Eq. 1.86 including all potential and boundary terms. One can prove 
that G(x,x',E) corresponds to the usual retarded Green function G r for < x, x' < a. 25 The 
continuum equivalents of Eqs. 1.82-1.85 are 

ip(a) = G(a,0,E)ihv L A, 

t = ihy/v^vE G(a, 0, E) = 2i^- Im Y, R G(a, 0, E) Im £ L , 
T = UmX R G r (a,0,E)lmX L G a (0,a,E). (1.89) 

I prefer linear algebra over differential equations and I don't have experience with the embed- 
ding formalism, so it won't be discussed further. It has been proven however that the formalism 
can be extended into a practical method for calculating transport through small systems, in- 
cluding all atomic details. 



1.5 Tight-binding 



In previous sections we have analyzed one-dimensional scattering problems starting from po- 
tentials and wave functions as functions of a continuous coordinate. We have discretized this 
representation to show how the scattering problem can be solved in a practical way. This ap- 
proach is natural if the potential variation is confined to the scattering region and outside this 
region the potential is constant. The incoming, reflected and transmitted waves are then simply 
plane waves. With the tunnel junction in mind, see Fig. 1.2, we would like to consider the 
situation in which the whole space is filled with atoms of one kind or another. In one dimension 
this gives an atomic wire, with atomic potentials everywhere along the wire, as is illustrated in 
Fig. 1.11. 

The usual approach is to construct a representation on a basis of atomic orbitals, i.e. expand 
the wave functions in fixed atomic orbitals Xi( x ) 



(1.90) 



where Xi denote the positions of the atoms. In chemistry this is known as the LCAO represen- 
tation (linear combination of atomic orbitals), in physics it is usually called the tight-binding 
representation. The wave function is represented by the column vector of the coefficients 



f \ ^ 

Ci-l 

Ci 
Cj+l 



V 



(1.91) 



/ 



Since we want to solve a scattering problem, the atomic orbitals should cover all of (one- 
dimensional) space, so the vector has an infinite dimension, i.e. i = — oo, ... , oo. Operators A 



25 G r (x, x', E) is defined by (^E + ir)- H^j G r (x, x' , E) = 6(x— x'), for — oo < x, x' < oo, with H the Hamiltonian 
without the boundary terms. 
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V(x) 




left region 




potential barrier I right region 



Figure 1.11: Top: atomic chain. The leads (left and right region) are periodic chains of iden- 
tical atoms. The middle region contains different atoms and/or disorder. Bottom: schematic 
representation of the potential along the chain. 



are represented by infinite dimensional matrices A with matrix elements Aij = (xi\A\xj) and 
the Schrodinger equation becomes 



(El — H) i/> = 0, 



(1.92) 



with H the Hamiltonian matrix. 26 To simplify the algebra we take just a single atomic orbital 
per atomic site. Moreover, we assume that the diagonal elements of H are Hii = hi] its off- 
diagonal elements are = /3j and all elements Hj^ = for j > i + 1. The other matrix 
elements are then set by demanding that the Hamiltonian matrix is Hermitian, i.e. = H. 



H 








V 
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PU 








i-1 


hi 


PI 








Pi 


hi+i 















\ 



(1.93) 



/ 



The model is called the nearest neighbor tight-binding approximation by physicists and among 
chemists it is known as the Hiickel approximation. The approximation is not essential and the 
formalism that I will explain below can be made to work for any LCAO representation, but here 
I want to keep the expressions as simple as possible. 27 

We divide our system into three parts: a left lead, a scattering region, and a right lead. The 
left and right leads are perfect materials. They consist of identical atoms at equal distances 
cll (for the left lead) and cir (for the right lead). In other words, the leads have translational 
symmetry. Matrix elements in the leads must be identical, i.e. hi = h L / R and /3j = Pl/r for the 
left/right leads. Only in the scattering region do we have site dependent matrix elements. The 
basic idea is illustrated in Fig. 1.12. 



26 I want to keep things as simple as possible, so I am using an orthogal basis. If the atomic orbital basis is 
non-orthogonal, then introduce an overlap matrix S, with matrix elements Sij = (Xi\Xj)i an d substitute I by S. 

27 In one-dimensional tight-binding one can assume that the hopping parameters fi i are real, without loss of 
generality. However, in three dimensions or in magnetic fields this need no longer be true. Therefore, to be 
prepared for the more general case, I keep complex fi^s for the moment. 
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Figure 1.12: Nearest neighbor tight-binding model of an atomic chain. The periodic left and 
right regions are characterized by the on-site and hopping matrix elements h L / R and The 
scattering region has site dependent matrix elements hi and 



1.5.1 Ideal lead modes 

Writing out Eq. 1.92 in its components gives 

-^_iCi_i + (E - fn)a - (3*Ci +1 = 0, (1.94) 

with i running from — oo to oo. These equations have the same mathematical form as the 
discretized Schrodinger equation of Eq. 1.58 if we make the substitutions 

h 2 h 2 

We are going to solve the scattering problem by mode matching. First we have to find the modes 
of the ideal leads. For sites i in the left lead the matrix elements are site independent and Eq. 
1.94 becomes 



-/9 L Ci_i + {E- h L )a - (3* L c i+1 = 0. (1.95) 

A little pondering shows that this equation is mathematically the same as a discretized Schrodinger 
equation for a particle in a constant potential. Since "the same equations have the same solu- 
tions" , the solutions must be (discretized) plane waves, i.e. 

c n = Ae ikLnaL , (1.96) 

where the distance between the atoms is the "discretization step". The same holds for the 
right lead, replacing the subscript L by R. 

One can give it a somewhat more mathematical flavor. The Bloch-Floquet theorem states 
that functions in consecutive cells of a periodic system are related by a constant amplitude/phase 
factor A, i.e. 28 

if Ci-i = c then Cj = Ac and q + i = A 2 c. (1-97) 



28 In the physics literature this is known as the "Bloch theorem", and in the mathematics literature as the 
"Floquet theorem" . Google a bit if you want to know the history. 
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Now is the time where we simplifying things a little bit by choosing (3 real, 29 so Eq. 1.95 becomes 



-0+(E- h)X - f3X 2 = 



(1.98) 



E-h 



E-h 



where you can add the subscript L/R for the left and right leads. The roots A can be given a 



more familiar form. For 



E-h 



2/3 



< 1 we define a wave number k by 30 
cos{kaj 



2(5 



which leads to the simple form 



(1.99) 



(1.100) 



A± is called the Bloch factor. Using Eq. 1.97 recursively, i.e. c n = X n co, then leads to Eq. 
1.96, with A = Co, as expected. It describes propagating waves, where A + describes a wave 
propagating to the right, and A_ a wave propagating to the left. 

> 1 one can define k by 



For 



E-h 



cosh(fta) 



E-h 



(1.101) 



and obtain 



if > l; 



(1.102) 



Both these cases describe states that decay either to the right or to the left. These are called 
evanescent states. They are not acceptable as solutions to the one-dimensional Schrodinger 
equation because one cannot normalize them (not even in the wave packet sense). However, we 
will have a use for them later on in three-dimensional problems. 



1.5.2 Mode matching 

We have the modes of the ideal leads, so we can match them to the scattering region, where the 
matrix elements hi and /3j in Eq. 1.94 are site dependent. We assume that the scattering region 
is localized in space, so i runs from 1 to N. The procedure we have used to solve the discretized 
Schrodinger problem in Sec. 1.4 can be copied with only some small modifications. Using Eq. 
1.100, Eqs. 1.61 and 1.64 read 

c_i = AX L \ + BX L ]_ = AX L \ + {c Q -A)X L ]_- (1.103) 

CN+2 = CN+lXR t+ . 

29 which is always possible in one-dimensional systems, provided one has spin degeneracy and time-inversion 
symmetry, which we have here. 

30 You may recognize E(k) = h + 2/3 cos ka as the dispersion relation describing an electronic band in the nearest 
neighbor tight-binding model. 
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The convention is to write Cj +n = A™q, where n is an integer (positive or negative), and let the 
± indicate waves propagating to the left and right, respectively. 31 These relations can be used 
to substitute Eq. 1.92 by 

(£I-H')V = q, (1.104) 

similar to Eq. 1.66. if) is a finite dimensional vector that contains the coefficients q in the 
scattering region plus those of the two boundaries, i.e. i = 0, . . . , N + 1. q is the "source" 
vector of length N + 2, whose coefficients are zero, except the first one 



P L A{\- L ]_-\- L ] + }. (1.105) 



H' is a finite (iV + 2) x (iV + 2) Hamiltonian matrix. All its matrix elements are identical to that 
of the original Hamiltonian matrix, Eq. 1.93, except for the first and the last diagonal element, 
which are modified to 

-^0,0 = h + T lL (E); 

h N+1 + E R (E), (1.106) 

with 

Z L (E) = P L X L ^_; 

E R (E) = (3 R \ R ,+. (1.107) 

These self-energies contain all the information concerning the coupling of the scattering region 
to the leads. As before, they are complex and energy dependent through Eqs. 1.99 and 1.100. 

We have substituted an infinite dimensional problem, Eq. 1.92, by a finite dimensional one, 
Eq. 1.104!! 32 Mathematically Eq. 1.104 is the same as Eq. 1.66. Again "the same equations 
have the same solutions", so according to Eq. 1.72, the transmission amplitude becomes 

W+1 (1.108) 




v L A 

with, as in Eq. 1.77, the velocities given by 

v L /R = -^^ImS L/fl (E). (1.109) 
In addition, the Green function expressions given in Sec. 1.4.1 remain valid. 



31 In the one-dimensional case, one always has A + = 1/A_. In the three-dimensional case, this relation does not 
necessarily hold. That's why it is important to keep separate track of the powers of A and the ± indices. 

32 For those of you who have some background in this, you might suspect that the technique we are using here 
has something to do with a technique known as partitioning. This suspicion is appropriate. Partitioning is 
usually applied to operators and, in particular, to Green functions. I am applying it to wave functions here. 



Chapter 2 

Electron transport in three 
dimensions 

In this chapter I want to upgrade the discussion of the previous chapter to problems in three di- 
mensions. The tunnel junction of Fig. 1.2, for instance, is clearly not a strictly one-dimensional 
device. In Sec. 2.1 we will generalize Landauer's formula for the conductance to three dimen- 
sions. I will do this in steps of increasing complexity, i.e. start with plane waves and square 
barriers, as in the one-dimensional (ID) case, and end with atomistic three-dimensional (3D) 
devices. The techniques that are introduced should allow us to study transport at the nanoscale 
in any small system. In a tunnel junction the dimension of the device across the junction (de- 
termined by the thickness of the insulator layer) is usually much smaller that the dimensions 
parallel to the metal insulator interfaces. The modeling then becomes easier if one assumes the 
latter dimensions to be infinite and one uses translation symmetry along the interfaces. The 
latter requires doing a bit of solid state physics. Another "device" is shown in Fig. 2.8 below. It 
consists of a mono-atomic wire sandwiched between two electrodes. Such wires have been made 
of gold atoms, for instance. Clearly the wire is finite in all three dimensions and the electrodes 
are infinite. Our modeling must be able to cope both with finite and with infinite dimensions. 

As in the previous chapter, Landauer's formula expresses the conductance in terms of the 
transmission of incoming electron waves. In three dimensions the transmission amplitudes con- 
stitute a matrix, the transmission matrix, whose calculation is our primary task. Practical 
algorithms are explained in Sec. 2.2. In particular, the mode matching techniques of Sec. 1.4 
are generalized to 3D. Since in 3D the derivations of the expressions become quite lengthy and 
technical, I won't give all details. They can be found in the literature. 

2.1 Landauer in 3D 

2.1.1 Conductance of model interfaces 

We start with a generalization of the square barrier potential of Fig. 1.3. A two-dimensional 
example is shown in Fig. 2.1. The potential V(x, y) is separable; the barrier is in the x-direction 
and the potential is independent of y. The straightforward extension to three dimensions is a 
potential that is independent of y and z. Such a potential landscape is a simple model for a thin 
layer of an insulator sandwiched between two metals, i.e. a tunnel junction. 
The Schrodinger equation 

h 2 

E^(r) + — V^(r) - V(r)tP(r) = (2.1) 
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x 

Figure 2.1: A square barrier in two dimensions. In the left and right regions the potential is 
constant, V(x, y) =V\. In the middle region the potential is also constant, V(x, y) = Vo, where 
Vq > V\. The incoming, reflected and transmitted waves are given by Ae ik r , J Be ik ' r and Fe ikr , 
see Fig. 2.2. 



is separable in Cartesian coordinates 

h 2 d 2 



E + 



2m dx 2 

and its solutions can be written as 



V(x) } + 



2m dy 2 2m dz 2 



i/)(x,y,z) = (j)(x)e ik y y e ik * z 



(x)e ik ir r ; k,| 





k 

k z 



y ' 



(2.2) 



where e* k H' r describes a free particle in the direction parallel to the barrier. The function cj){x) 
is the solution of the one-dimensional scattering problem 



Ej. 



2m dx 2 



V(x) 



(x) = 0, 



(2.3) 



with the energy 



21.2 



E^r — E 



2m 



(2.4) 



Note that kn is a good quantum number 1 , which together with the energy E fixes the wave 
function. We say that kn defines a mode of the system. 

A view of the scattering geometry in the (A^kii) plane is given in Fig. 2.2. In the left and 
right regions the wave function is 



V> k|l ( r ) 



A 



e 2k-r _|_ r , (E)e lW ' r ; r in left region 



= ik'r 



.4 



tUE)e 



ikr 



; r in right region 



(2.5) 



1 a, set of two quantum numbers, actually. 
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Figure 2.2: Scattering from a planar square barrier as viewed in the (A^kii) plane. 

with k = (fca;,kj|) and k' = (— k x ,k^). I leave you the job of finding the expressions for the 
transmission and reflection amplitudes t^(E) and r^(E). In three dimensions the probability 
current is a vector 

J(r, t) = ^ [$(r, t) W*(r, t) - **(r, t)W(r, *)] (2.6) 

As in the one-dimensional case, for a stationary problem, \I/(r, t) = ^(r)e _ R" B *, the current is 
constant. In three dimensions it is actually more appropriate to use the phrase "current density" 
for J (as in classical electrodynamics), but I keep on using the phrase "current" for short. 



The key point 

For the experiments we are interested in, only J x matters, i.e. the x-component of the current. 
For devices like that shown in Fig. 2.1 (and all other devices that we will consider), one attaches 
macroscopically large electrodes to the left and right regions and applies a bias voltage between 
those electrodes. Only the current along the x-direction is then measured. 2 ' 3 So we are interested 



m 



'J:r 



2m 



dj>*{r) 
dx 



dx 



(2.7) 



From Eq. 2.5 one can easily show that the transmitted current carried by one mode ky is given 
by 



Jta 



*k,| (E) 



(2.« 



with the density p = \A\ 2 and the velocity in the x-direction v x = ^ 3L . Deriving an expression 
for the conductance follows the same steps as in Sec. 1.2. Applying a small bias AV = —ell 



2 For those of you who are experts in scattering phenomena, note that this is different from what you are used 
to in angle resolved three-dimensional scattering experiments in free space. There one is usually interested in all 
three components of the current. In addition, if the scatterer is localized in space, one applies a transformation 
to spherical coordinates. 

3 More complicated measurements are possible. In "multiprobe" experiments more than two electrodes are 
attached to the device. This is often done in combination with applying an external magnetic field, as in measuring 
the Hall effect. The Landauer formalism can be extended to include multiprobe measurements. This extension is 
due to Biittiker. 
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between left and right regions, the transmitted current is carried by all modes that have an 
energy E^ in the range from Ep — AV to Ep, see Fig. 1.7. 



J T = 2p 



E F -AV<E k <E F 



(2.9) 



where the factor 2 accounts for the degeneracy. We can use the same tricks as in Sec. 1.2. 
Using states that are normalized in a 3D, we have p = ^p-; compare to Eq. 1.32. 4 We write 
IjE F -Ay<E k <E F = Sfcz' wri ere one has to sum only over those states that have their 
energy in the indicated interval. Convert Y2k a one-dimensional integral, use v x = ^^rS 

and assume that (E) is independent of the energy in the small energy range AV. The algebra 
is 



Ji 



*-|| k x J 

L 2 ^ eKi(^ 



1 1 
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1 AV 



2 1 d£L 



l l ^ r E F-^t 

h dk x """ x L 2 irh f-^ Jp^^av^^I 

Km ^ 2m 



2 1,2 



tki, (-Ea; + 



2m 



dE x 



(2.10) 



The transmitted current is It = J Jt • dS = JpL 2 and the conductance Q = It/U is then given 
by 



^ = ^Ek(^) 2 = ^w 



(2.11) 



The expression is the Landauer formula, see Eq. 1.29, with the total transmission T expressed 
as a sum over the transmissions of the individual modes ky. One has to sum over the ky that 
contribute to the transmission at the Fermi energy Ep. 



The Fermi surface 

The Fermi surface is defined by the relation E^ = Ep, the Fermi energy being a materials 
constant. It can be visualized as a surface in reciprocal space, i.e. in three-dimensional k-space. 
For free electrons or electrons in a constant potential 




the Fermi surface is the surface of a sphere, as shown in Fig. 2.3(a). The Fermi surface can 
help to visualize which modes contribute to the transmission. The latter can be enumerated by 
projecting the part of Fermi surface with k x > (a hemisphere in this case) onto the ky = (k y ,k z ) 
plane, as shown in Fig. 2.3(b). All the modes ky within this projection exist at the Fermi energy 
Ep and they contribute to the transmission in Eq. 2.11. The scattering geometry in real space 
can be deduced from Fig. 2.2. ky = (0, 0) 5 corresponds to a wave with normal incidence to the 

4 The function ■^■e lk H ' r is normalized in a 2D box of size L, as you can easily check yourself. The function 
4>{x) is treated as in Eq. 1.32 and gets the normalization factor j^. The product gives the normalization factor 
of ip(x, y, z), see Eq. 2.2. 

called the T-point in solid state physics folklore. 
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x (a) 



(b) 



(c) 



Figure 2.3: (a) The Fermi surface, as defined by £k = Ep, of electrons in a constant potential 
is a sphere, (b) The projection of the surface in the kn = (k y ,k z ) plane. The shaded area 

denotes all the ky modes that contribute to the transmission at Ep. (c) The transmission t^.. 
as function of kn. Red indicates a high transmission and blue a low transmission. The highest 
transmission is for kn = (0,0). It decreases to towards the edge of the circle. 



barrier. The larger fen = knl, the more glancing the incidence of the corresponding wave. At 

Ep, so k x = 0, and the corresponding wave 



h h 2 

2mH| 



the edge of the circle in Fig. 2.3(b) one has 
propagates parallel to the barrier. 

For a simple square barrier one can easily calculate the transmission. If the Fermi energy is 
less than the barrier height, the transmission of a single mode is given by 



2l.2 



iky (Ef) = Tid{E x ) = Tip>{Ep 



2m 



(2.12) 



with T\p, given by Eqs. 1.47 and 1.42. The result is visualized in Fig. 2.3(c). Since T\p> is a 
monotonically increasing function of the energy, the maximal transmission is for kn = (0, 0), i.e. 
for normal incidence. The transmission decreases monotonically with increasing fen, i.e. with 
the angle of incidence, until it is zero for parallel incidence. Such a simple Fermi surface and a 
simple transmission are typical for free electrons. For real materials both the Fermi surface and 
the transmission are much more complex. I will come back to in Sec. 2.1.4. 



2.1.2 Conductance of model wires: the ballistic conductance 

The way we have derived Eq. 2.11 means it can easily be generalized to any system whose 
Hamiltonian is separable in an x-term and a yz-term. For instance, the wave functions of the 
wire shown in Fig. 2.4 can be written as 

i/)(x, r±) = cf)(x)r] n (r ± ), (2.13) 

where n is a set of two quantum numbers labeling the modes; see Eq. 2.2. The conductance of 
this quantum wire can be expressed as the sum of the transmissions of the individual modes 




(2.14) 



n 



Suppose we are dealing with an ideal wire in which all the modes are fully transmitted, 
which is the definition of an ideal wire. 
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Figure 2.4: A simple (quantum) wire. 



The transmission is maximal, i.e. \t n (EF)[ 

Gba\ 



1 for each mode n, and the conductance is 

2 M(E F ), 



e 



(2.15) 



where M(Ep) is the number of modes at the Fermi energy, supported by the wire. £?bai is called 
the ballistic conductance. The number of modes at fixed energy in a wire of finite cross 
section is finite. So even for a perfect wire the conductance is finite!! If the cross section is 
macroscopically large, the number of modes is very large. Measuring the ballistic conductance 
experimentally is then impossible. A real wire always contains impurities and imperfections (lat- 
tice defects, grain boundaries, impurity atoms, etc.), and scattering at these defects dominates 
the conductance. However, thin and small wires can be made without any defects. Note that, 
since the number of modes is an integer, the conductance of Eq. 2.15 is quantized. 

A very clear example of quantized conductance is presented by the experiment of van Wees 
et al. on the quantum conductance of a channel in a two-dimensional electron gas (2DEG). A 
2DEG is formed at an interface between two well-chosen semiconductors. Putting electrodes on 
top of the 2DEG it is possible do define a one-dimensional channel (the "wire") by a suitable 
electrostatic potential profile, as illustrated in Fig. 2.5. The width of the channel determines 
the number of modes at the Fermi energy M(Ep). Upon widening the channel, which is done by 
making the gate less repulsive for electrons, the number of modes increases and the conductance 

2 _ 

increases. However, it increases in steps of This is demonstrated by the experimental results, 
shown in Fig. 2.5. 

2.1.3 Conductance of combined interfaces and wires 

In the real world interfaces nor wires are infinite. This leads to Hamiltonians that are not 
separable in an x-term and a yz-term. A two-dimensional example is shown in Fig. 2.6, which 
represents a canyon through a 2D square barrier. It is a simple model for a finite wire connected 
to two electrodes (the left and right regions). An incoming wave Ae l ^' r with k = {k x , kn) can be 
scattered into any other wave Ce* k ' r with k' = (k' x , k[,) provided the energy is conserved, i.e. 



h\kl + kl) 



E 



2m 2m 
The idea is shown in Fig. 2.7. The wave function of Eq. 2.5 is generalized to 



(2.16) 



VUr) 
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+ £ k '.' r k ll .k"( £; ) e 



ik r 



A 



£k' '„ V (/-> 



ik' 



; r in left region 



; r in right region 



(2.17) 



where k' labels the transmitted waves, i.e. k' x > 0, and k" labels the reflected waves, i.e. 

all at the same energy. The reflection and transmission amplitudes r, , » and t, , / indicate 

K||,K|| K||,K|| 

the probability amplitudes of scattering from a mode kn to a mode k or k . 
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Figure 2.5: Quantized conductance of a ballistic waveguide Left: top view of the device. The 
current flows from the source to the drain electrode. The gate has a negative potential and 
repels the electrons towards the middle. This leaves an effective channel for the electrons. At 
a small gate potential the channel is wide (indicated in red), and at a large gate potential the 
channel is narrow (indicated in blue). Right: measured conductance as function of the gate 
potential. The conductance is quantized in units of e 2 /irh, as shown by the plateau's; M is the 
number of modes at Ep contributing to the conductance. See: B. J. van Wees et al, Phys. Rev. 
Lett. 60, 848 (1988). 



The derivation of the Landauer formula follows the same steps as in the previous section. 
For instance, the current in the x-direction, carried by ip^ (r) is determined from Eq. 2.7 



'kii.k 



(E) 



(2.18) 



with p = \A\ Z and v' x = compare to Eq. 2.8. The generalization of Eq. 2.9 is 



E F -AV<E k <E F k ' 



'kn,k 
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'k,,,k 



(2.19) 



where the sum over the energies has been replaced by a sum over the states that have their 
energy in the right interval. Note that one does not have to sum over k x . By fixing k[|,k,i and 
k' x one has fixed k x , because of Eq. 2.16. One can follow the steps of Eq. 2.10 and obtain an 
expression for the conductance which generalizes Eq. 2.11 



g=-y 



k ,1c. 



kii,k 



{E t 



(2.20) 



The transmission amplitudes k / are collected in a transmission matrix t, which allows the 
conductance to be expressed in a compact form as 



£ = ^Tr 
nh 



nh 



T{E F ). 



(2.21) 



This is the Landauer formula for non-uniform wires and interfaces. 
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Figure 2.6: A simple model for a wire connected to two electrodes (the left and right regions) 
as a canyon through a square barrier. The incoming wave Ae ik ' v can be reflected to any wave 
i?e ik r or transmitted to any wave Fe 1 ^ r , provided the energy is conserved. 




Figure 2.7: Scattering from a finite wire between two electrodes as viewed in the (k x , kii) plane. 
Only one of the possible reflected (k", kj|^ waves and one of the transmitted (k' xl kj^ waves are 
indicated. 

2.1.4 Conductance of atomistic wires and interfaces 

We derived the Landauer formula, Eq. 2.21, using a model in which the incoming, reflected 
and transmitted waves in the asymtotic regions are plane waves e lkr . 7 One might wonder how 
much of this derivation depends on having plane waves. Not too much, actually. We need to 
know what the waves in the asymtotic region look like, but they do not need to be plane waves. 
Consider the example shown in Fig. 2.8. It shows a mono-atomic wire of three atoms between 
two planar electrodes. Such wires have been made of gold and other metals (albeit with a less 
ideal structure). 8 The scattering region comprises the wire and a region around the wire in 
which the atomic potentials are different from bulk atomic potentials. Exactly how large that 
region is, depends somewhat on the system, but it will be clear that far into the electrodes, i.e. 
far into the left or right leads, the atoms are bulk metal atoms. 9 



7 The regions outside the scattering region, i.e. the left and right regions are called the asymtotic regions. 
8 by a technique called "mechanical break junctions", for instance, see J. van Ruitenbeek. 

9 For metallic electrodes the potential at a few atomic layers from the surface is usually indistinguishable from 
the bulk potential, because screening effects in metals are large. 
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incoming wave 




Figure 2.8: Schematic mono-atomic wire consisting of three atoms between two metal electrodes. 
Outside the scattering region one considers the electrodes as bulk metals, whose electronic states 
are represented by Bloch waves. 

Bloch modes 

Assuming that the bulk leads consist of perfect crystalline material, the Bloch-Floquet 
theorem tells us that the electronic states are represented by Bloch waves, which have the 
form 

0kn(r) = ^kn(r)e^ r . (2.22) 

Here u^ n (r) is a periodic function in three dimensions, and the primitive unit corresponding to 
the three periods is called the unit cell, k is a wave vector within the reciprocal unit cell (called 
the first Brillouin zone) and n is the band index (a quantum number that distinguishes between 
states with the same wave vector). Proof of this can be found in any book on solid state physics. 
These Bloch waves replace the plane waves as modes in the scattering formalism. 

The relation between energy and wave vector E^ n is called a dispersion relation. It can 
be obtained by solving the Schrodinger equation for electrons in the periodic crystal potential. 
As before, the Fermi surface is defined by the relation E^ n = Ep, which gives a surface in 
reciprocal space, i.e. in k-space, as in 2.3(a). Actually, it gives a surface for each n, called a 
sheet of the Fermi surface. 

Fig. 2.9 gives a few examples of Fermi surfaces for different metals. The Fermi surface of 
some metals resembles that of free electrons; compare Figs. 2.9(a) and 2.3(a). These are called 
"simple metals" ; they usually involve s- and/or p-electrons only. In other cases the Fermi surface 
can be quite complicated. It can consist of more than one sheet and has a shape that is far from 
free electron like (i.e. spherical), as is shown in 2.3(b). This is common among metals involving 
d-electrons. 

The Landauer formula is still valid, however. The conservation of energy, Eq. 2.16, has to 
be replaced by the more general relation 

E\ sn = E = E\sj n i, (2.23) 

which basically states that one can scatter from a state kn to any state k'n', provided both points 
are on the same energy surface (the Fermi surface). In Eq. 2.17 one has to replace the plane 
waves e lkr ,e lk r ,e lk ' r by Bloch waves 0kn( r )> 0k'n'( r )> 0k"n"( r ) an< ^ the scattering amplitudes 
now also contain the band index, i.e. r, , » „ and L . / ,. Another vital ingredient needed 

' k||Tj,k|, n" k||n,k,|n' ° 
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Figure 2.9: (a) The Fermi surface of Cu (copper). The wireframe indicates the first Brillouin 
zone (the reciprocal unit cell). The Fermi surface is almost spherical with a few holes punched 
through it. In Cu there is only one s-like band crossing the Fermi energy, which gives rise to 
one sheet, (b) The Fermi surface of fee Co (Cobalt) for the minority spin electrons. Co is a 
ferromagnetic material. Only one s-like majority spin band crosses the Fermi energy, which 
gives rise to a Fermi surface for majority spin electrons that resembles that of Cu. Three d-like 
minority spin bands cross the Fermi energy, giving rise to three sheets of the Fermi surface that 
are far from spherical, as is shown here. 



is an expression for the velocity v' x . In solid state textbooks it is shown that the Bloch velocity 
is given by 



1 dEv 



h dk' 



(2.24) 



A little pondering shows that this is just the relation we need to complete the derivation of the 
Landauer formula. 
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(2.25) 



An example: the conductance of an interface between copper and cobalt 

Does it matter whether one takes real metals with Bloch modes as electrodes, or jellium metals, 
where the modes are plane waves? As an example that it does matter I will discuss the con- 
ductance of a simple interface between the metals copper (Cu) and cobalt (Co). Bulk Cu has 
a fee (face centered cubic) close packed structure. Its (111) surface is close packed and it forms 
the crystal termination with the lowest surface energy. Bulk Co has the hep (hexagonally close 
packed) structure. However if one deposits a thin layer of Co on Cu(lll), then Co assumes a 
fee structure with the lattice parameter of Cu. So the (111) interface between Cu and Co has 
a simple structure. So far we have considered spin degenerate systems, i.e. spin-up electrons 
are described by the same wave functions as spin-down electrons. This hold for Cu, which is a 
non-magnetic material. It does not hold for Co, which is a magnetic material, and both spins 
have different wave functions. In addition the number of spin- up electrons is not equal to the 
number of spin-down electrons (which is why the material is magnetic). The obvious names 
majority spin and minority spin are used. 
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The Landauer formula can be generalized to include the spin states explicitly in a straight- 
forward manner 



e 2 



2-kH 

k|in<7,k|| n'a' 



^kii7i<7,k{| n'a' ' / 



2 

(2.26) 



where a, a' = — 2 - ^ n & bsence of spin-orbit coupling the scattering potential does not flip the 
spin and the transmission matrix is diagonal in the spin t, , > , , = t, , / , b nn i . This is 

1 ° 1 k||ncr,k||n'cr' ky ncr, ky n'tr 00 

accurate for (magnetic) 3d transition metals and lighter elements. 

One can use the same procedure as in Fig. 2.3 to visualize the transmission. The Fermi 
surfaces of Cu and of the minority spin of Co are given in Fig. 2.9. The Fermi surface of the 
majority spin of Co resembles that of Cu. One has to project these surfaces onto the ky plane, 
i.e. along the transport direction, which is the (111) direction in this case. In Fig. 2.9(a) (111) 
is the direction pointing from the origin to the center of a hexagon. The projection of the Fermi 
surface of Cu is shown in Fig. 2.10(a). The projected Fermi surface for the majority spin of Co 
is shown in Fig. 2.10(b). As stated, it is quite similar to that of Cu. The Fermi surface for the 
minority spin of Co is much more complicated, as shown in Fig. 2.9 (b). Projecting all three 
sheets in the (111) direction on top of one another gives the result shown in Fig. 2.10(c). 

For a perfect crystalline interface ky is a good quantum number and the transmission matrix 
is diagonal in ky (as well as in a) * k k ' n , a , = t^na^n'a^^' S aa ' . For the majority spin 

2 

tu.,, io-,k|, la describes the transmission from the state kyla on the projected Fermi surface of 
Cu (indicated by a black dot in Fig. 2.10(a)) to a state kiilcr on the projected majority spin 
Fermi surface of Co (indicated by a black dot in Fig. 2.10(b)). 10 The calculated transmission 

2 

tu.,, i<7,k|i i<t is shown in Fig. 2.10(d) in the ky plane. The transmission is a number between 
and 1. Obviously the state knltr has to exist on the Cu Fermi surface; otherwise one has no 
electrons to start with. Around ky = (0, 0) the Cu Fermi surface has a "hole" where there are 
no states, which is indicated in white in 2.10(a) and (d). Then the state ky la has to exist on the 
Co Fermi surface; otherwise there is no state to transport the electrons to and the transmission 
is 0. These are the blue area's in Fig. 2.10(d). The remaining area is almost entirely red, which 
means that there the transmission is close to 1. 



Fig. 2.10(e) shows the calculated transmission Yl n '=i *k||i<T,k|in'<T between the minority spin 
modes of Cu and Co. As one can observe, the transmission pattern is much more complicated 
than in the majority spin case with a fine-grained variation of the transmission between and 
1. Now compare Figs. 2.10(d) and (e) to Fig. 2.3 (c). It represents the difference between using 
real metals and using jellium. 11 

Integrating over the area shown in Fig. 2.10(d) gives the majority spin conductance; the 

2 

calculated number is £? ma j = ^0.73. Integrating over the area shown in Fig. 2.10(e) gives 

2 

the minority spin conductance Gmin = 2^0-66. The difference in the conductance between the 
majority and the minority spin modes is 10%. One can enlarge this difference by making a 
multilayer of alternating Cu and Co layers since such a multilayer then contains many Cu/Co 
interfaces. The majority /minority difference in the conductance essentially gives the so-called 
GMR (giant magneto-resistance) effect that is measured in such multilayers. 12 



10 n = 1 and n' = 1, since both Fermi surfaces have only one sheet. 

11 Jellium also does not give magnetism at these electronic densities, but I don't want to make a list of what's 
wrong with jellium. 

12 The GMR effect is used to make sensitive sensors for magnetic fields, which are used in the heads of magnetic 
hard drives. It contributed to the rapid development in hard disks to the >10 2 Gb disks we have today. 
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Figure 2.10: (a) Fermi surface of Cu projected in the (111) direction on the k| plane, (b) Fermi 
surface of the majority spin of Co projected on the same plane, (c) The three sheets of the Fermi 
surface of the minority spin of Co projected in the same plane. Blue, green and red indicate 
area's where 1,2 and 3 states are projected onto the same kii point, (d) The transmission across 
a (111) interface between the Cu and the majority spin in Co, as a function of kn. The color 
scale indicates the transmission from (blue) to 1 (red), (e) As (d), but for the transmission 
between Cu and the minority spin of Co. The white area's indicate where there are no states in 
Cu. 



2.2 Calculating the transmission 



As will be clear by now, the key quantity for calculating the conductance is the transmission 
matrix t. This is quite a bit more complicated than in the one-dimensional case, but in principle 
all the techniques discussed in Sees. 1.4 and 1.5 can be extended to calculate transmission 
matrices in three dimensions for systems containing real atoms. Since the tight-binding example 
given in Sec. 1.5 is chemically the most intuitive, I use this and focus upon the mode matching 
technique. 

The principal idea is to divide your system into layers of atoms, normal to the transport 
direction. One chooses the layers to be sufficiently thick such that the Hamiltonian matrix only 
contains matrix elements that either couple atoms within one layer, or atoms that are in nearest 
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Figure 2.11: Hamiltonian of a tunnel junction divided into layers. The transport direction is 
along the horizontal. The left (L) and right (R) leads are ideal periodic wires containing the 
layers i = — oo, . . . , and i = S + 1, . . . , oo, respectively. The layers i = 1, . . . , S consitute the 
scattering region. 



neighbor layers. 13 The Hamiltonian matrix of Eq. 1.93 then becomes 
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(2.27) 
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A schematic representation of the structure of the Hamiltonian is given Fig. 2.11. The matrices 
Hj contain the interactions between atoms within the layer i, with i = — oo,. . . , oo, whereas 
the matrices Bj describe the coupling between the layers i and i + 1. Both of these are N x N 
matrices, where N is the total number of atomic orbitals for all atoms in a layer. The scattering 
region is localized in the layers % = 1, . . . , 5. 14 

This representation is valid both for systems with a finite cross section (i.e. wires), and for 
infinite layered systems that are periodic along the interfaces (i.e. along the layers). In the last 
case N refers to the number of atoms within the unit cell, kii is then a good quantum number 
and the matrices Hj(kii) depend on this quantum number. I won't discuss the exact expressions; 
they can be found in the literature. In fact, in order to simplify the notation I will omit the 
quantum number kii from now on. I will use the phrase "wire" to indicate both systems with a 
finite cross section and systems with infinite periodic interfaces. 

The wave function of Eq. 1.91 is generalized to 



/ 



\ 



(2.28) 



13 The formalism can be extended to include interactions of a longer range, but the expressions become a bit 
messy. So far there has been no need for this extension. 

14 Note that compared to Sees. 1.4 and 1.5 I have a slight change in notation here. Now is the dimension of 
the basis within one layer, and S is the number of layers in the scattering region. Sorry about that; it's what you 
get if you merge texts that have a different origin. 
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where Cj is the N- vector of coefficients of the atomic orbitals of all atoms in a layer. As in Sec. 
1.5 we divide our system into three parts, with i = — oo, ... ,0 corresponding to the left lead 
(L), i = 1, . . . , S to the scattering region (S) and i = S + 1, . . . , oo to the right lead (R). Using 
the mode matching technique the scattering problem is solved in two steps. In the first step the 
Bloch modes of the leads are calculated and in the second step these are then matched to the 
scattering region. 



2.2.1 Ideal lead modes 

The leads are assumed to be ideal wires characterized by a periodic potential. It is then ap- 
propriate to identify a layer with a translational period along the wire. By construction, the 
Hamiltonian matrix then is the same for each layer in the leads, i.e. Hj = H^/^ and Bj = Hl/r 
for the left/right leads, see Fig. 2.11. Eq. 1.95 is generalized to 

-B LARCi _i + (El - H L/R ) Ci - B[ /R c l+1 = 0, (2.29) 

where 1 is the N x N identity matrix. 15 We make the same ansatz as in Eq. 1.97, namely that 
the coefficients in successive layers are connected by a Bloch factor A 

Cj_i = c; Cj = Ac; c i+1 = Acj = A 2 c (2.30) 

to get the equation 

—Be + X(E1 — H)c — A 2 B^c = 0, (2.31) 

where the subscripts L/R have been omitted to simplify the notation. Remember that we work 
at a fixed energy E, so Eq. 2.31 is a quadratic eigenvalue equation in A of dimension N. There 
are standard tricks to solve such equations. For instance, by defining d =Ac one can convert it 
into 



— B £71 — H 7 Bt 



1 \ . ( 1 

A 



C d ) =0. (2.32) 



This is a linear (generalized) eigenvalue problem of dimension 2N, which can be solved using 
standard numerical routines. 

It can be shown that this equation indeed generally has 2N solutions, which can be divided 
into N right-going modes and N left-going modes, labelled by a "+" and a "— " subscript as 
in Eqs. 1.100 and 1.102. Right-going modes are either evanescent waves that are decaying to 
the right, or waves of constant amplitude that are propagating to the right, whereas left-going 
modes are decaying or propagating to the left. Figs. 2.12 and 2.13 give you a simple idea. 

In contrast to the ID case, we find in three dimensions at a fixed energy in general both 
evanescent and propagating modes. We denote the eigenvalues and eigenvectors of Eq. 
2.31 by 16 

A±,„; u±, n ; n = 1, . . . , N. (2.33) 

Together these states form a complete basis set. In the following we assume that the vectors 
u± n are normalized. However, note that in general they are not orthogonal. 17 One can easily 



15 Again, for non-orthogonal basis sets one can substitute 1 by S, the overlap matrix. 

16 Again I remind you that I have omitted the ky quantum number. So if you see the index n to denote a mode 
in the following, you can substitute it by kyn or by kyncr to denote the full mode. 

17 This follows from the quadratic eigenvalue problem, Eq. 2.31, or the equivalent generalized linear one, Eq. 
2.32. 
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Figure 2.12: Right-going, i.e. "+", modes. Black curve: propagating mode, |A+| = 1. Blue and 
red curves: examples of evanescent modes, |A+| < 1. 




Figure 2.13: Left-going, i.e. "— ", modes. Black curve: propagating mode, |A_| = 1. Blue and 
red curves: examples of evanescent modes, |A_| > 1. 

distinguish right- from left-going evanescent modes on the basis of their eigenvalues. Right-going 
evanescent modes, which decay to the right, have |A +jn | < 1 and left-going evanescent modes, 
which decay to the right, have |A_, n | > 1, see Eq. 2.30 and Figs. 2.12 and 2.13. Propagating 
modes per definition have |A± n | = 1, so here one has to determine the Bloch velocity in the 
propagation direction and use its sign to distinguish right from left propagation. One can show 
that for a tight-binding form of the Hamiltonian, the general expression for the Bloch velocities 
of Eq. 2.24 becomes 



2a T 

v±, n = — r Im 
n 



A±, n i4 Btu±J , (2.34) 



where a is the thickness of the layer. The derivation of this expression is rather technical and 
can be found in the literature. In addition one can show that the Bloch velocity is non-zero only 
for propagating modes, i.e. propagating states have a Bloch velocity equal to zero. 18 

Since the eigenvectors are non-orthogonal, it is convenient to define dual vectors u±, n by 

U± n U ±,m = S n ,m', U± n™±,m = #n,m- (2.35) 



18 Which is what you might expect. Evanescent modes contribute a zero current, i.e. no particles, energy, or 
whatever is transported by evanescent modes. 
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Any wave function in the leads can be expressed as a linear combination of the lead modes. This 
can be done in a very compact way by defining the two N x N Bloch matrices for right- and 
left-going modes 

N 

F± = ^A^u^G^. (2.36) 

n=l 

These are the generalizations of the Bloch factors of the one-dimensional case, see Eqs. 1.100 
and 1.102. Note that one can easily construct powers of Bloch matrices 

N 

n = E A ±,- u ±>U- ( 2 - 37 ) 

n=l 

This is valid for any integer i, due to Eq. 2.35. A general solution in the leads can now be 
expressed as a recursion relation 

d = c +;l + c_ ;l = F^'c+j + Fr j c_j. (2.38) 

In a scattering problem one usually fixes the coefficients in one layer using boundary conditions. 
By using Eq. 2.38 one can then determine the solution in all the layers of the leads. 

2.2.2 Mode matching 

One can use these recursion relations to set up a set of equations that properly match the leads 
to the scattering region. The scattering region is defined by the layers i = 1, . . . ,S, see Fig. 
2.11. Immediately left of the scattering regions one has the recursion relation 

c-i = F-^c+,0 + F2 ) 1 _c_, . 

Writing c_ ; o = Co — c +i o one gets 

c+,o + F^_c . (2.39) 

The scattering boundary condition is introduced in the following way. The vector c +i o is treated 
as the source, i.e. as the incoming wave from the left lead, e.g. a specific (propagating) mode 
of the left lead. 

c+,o = UL !+ , m . (2.40) 

Immediately right of the scattering region one has 

cs+2 = Ffl )+ c +) s+i + F R _c_ s + i. 

Here the scattering boundary condition states that we have no incoming wave from the right 
lead, i.e. c_ ; s + i = 0, so 

cs+2 = Ffl.+cs+i. (2.41) 
The relations Eqs. 2.39 and 2.41 can now be used to "simplify" the tight-binding equations 

-B iCi _i + (El - Hi)a - Bjc i+ i = 0. (2.42) 



c-i 



-F 



L,- 
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Ordinarily i runs from — oo to oo, but for i = we can write 
-B L c_i + (El - H L )c - Blci = 0^ 



B, 



F7 -FT 



(El - H L - BlF^Jco - B[ci 
Likewise one can write for % = S + 1 

B 5 c 5 + (£71 - H i? )c s+ i - B^cs+2 
-B s c s + (El -H R - B^F i?i+ )c 5 +i 

Eq. 2.42, for i = 1, . . . ,5, and Eqs. 2.43 and 2.44 can be collected into 

(£1-H')V = Ql,+,™, 



(2.43) 



0^ 
0. 



(2.44) 



(2.45) 



with 



/H L + E L (E) B\ 



H' 



B: 








B 



-i 



Bj_i Hj B^ 













Bt 



B, 

B S H R + t R (E) ) 



(2.46) 



and 



/ c \ 



/b, 



; Qi,+,r. 



V c s+i J 







UL,+,m 



(2.47) 



The quantities 



^l(E) 
Xr(E) 



BlF^ 



(2.48) 



are called the self-energies of the left and right leads, respectively, just as in Eq. 1.107. They 
contain all the information about the coupling of the scattering region to the leads, as well as the 
information about the scattering boundary conditions. Note that the self-energies depend upon 
the energy E, since are expressed in the Bloch matrices and thus in the lead modes and the 
latter have been determined at a fixed energy E. Moreover, the self-energies are non-Hermitian, 
which makes the Hamiltonian H' non-Hermitian. 19 An important observation is that Eq. 2.45 
has a finite dimension! 20 

Eq. 2.45 represents a set of linear equations of dimension (5 + 2) x N. These can be solved 
using the common techniques, e.g. Gaussian elimination (or LU decomposition, as it is also 



19 This is logical if you know your quantum mechanics. A particle in the scattering region has a finite probability 
to leak into the leads, i.e. to disappear from the scattering region. In other words, the particle has a finite lifetime. 
This can only be expressed by non-Hermitian "Hamiltonians" , as Hermitian ones would only give real energies 
and time factors e~T Bt yielding constant probabilities, i.e. infinite lifetimes. 

20 If you suspect that the "partitioning technique" is behind this, you are right. 
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known), followed by back substitution. Since the matrix H' has a block tridiagonal form (the 
blocks being of dimension N), one can make use of this special form to make the Gaussian 
elimination algorithm efficient. The details can be found in the literature. Transmission 
matrix elements are obtained by expanding the wave function in the right lead into modes 21 

N 

C S+1 = U R,+,nt'n,m => 4,m = Ufl,+, n C S+l ( 2 - 49 ) 
n-1 

and normalize them with the velocities to ensure a unitary scattering matrix, see Eqs. 1.53 and 
1.72. 

tn ' m = ^^77 5 k+,n c S+i (2.50) 

The velocities are calculated from Eq. 2.34. 

By letting m in Eqs. 2.40 and 2.47 run over all (propagating) modes u^ !+ m of the left lead, 
all transmission matrix elements can be obtained. This sounds actually more complicated than 
it is. Gaussian elimination (LU decomposition) is the most time consuming algorithmic step 
in solving the linear equations of Eq. 2.45. This step however is independent of QL,+,m- The 
latter only enters in the second algorithmic step, i.e. back substitution, which is much cheaper 
and is performed for each QL,+,m- Moreover, the back substitution step does not have to be 
complete, since only the coefficients cg + i are needed to calculate the transmission. 22 

Although the mode matching formalism allows for calculating t' n m matrix elements according 
to Eq. 2.49 between all states n and m, i.e. propagating states and evanescent states, only 
propagating states contribute to the transmission. In this sense it is allowed to restrict n and 
m to the propagating states only. This does not mean that one can throw out the evanescent 
states altogether, though. One needs to include them in properly matching the scattering region 
to the leads, cf. Eq. 2.48. In other words, the Bloch matrices, Eq. 2.37, need to be constructed 
from all N modes, propagating and evanescent. 23 Both of them are needed for a complete basis. 



2.2.3 Green function expressions 

For those of you who cannot live without Green functions, I can give you the Landauer formalism 
expressed in Green functions. The derivations are rather technical, so I will skip most of them. 
You can find them in the literature. With respect to the Hamiltonian matrix of Eq. 2.46 a finite 
Green function matrix can be defined as 

G(E) = (El — H') -1 . (2.51) 

It can be calculated by matrix inversion using essentially the same block Gaussian elimination 
scheme as in solving the set of linear equations, Eq. 2.45. The matrix H' is non-Hermitian; its 
eigenvalues are not real, so the Green function matrix can be evaluated for real energies. 24 

With respect to the original infinite Hamiltonian H of Eq. 2.27 one can define the usual 
retarded (infinite) Green function matrix 

G r (E) = [(E + i V )l~-H}- 1 , (2.52) 

21 Again, replace n by kyn or kiincr, when appropriate. 

22 However, if one wants to have the full wave function ip or calculate the reflection coefficients, the back 
transformation has to be complete. It is not expensive, though (compared to the Gaussian elimination step). 
23 0nly the full set forms a complete basis. 

24 The exceptions are truly bound states of H'. These occur at real energies, at which the Green function matrix 
thus has poles. However, these occur outside the energy range that is of interest for scattering. At bound state 
energies only evanescent modes exist in the leads, so this situation is easily identified in practice. 
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where one needs the infinitesimal 77 to avoid the poles on the real axis. The advanced Green 
function G a (E) can be obtained from Eq. 1.80. One can show that for layers in the scattering 
region one has 

Gij(E) = G r id (E); i,j = 0,...,S + l. (2.53) 

if one properly takes the lim^o- in terms of the Green function, the wave function in the 
scattering region can be written as 

i> = G(E) Q L>+ , m , (2.54) 

and from Eq. 2.49 one obtains 



tn,m = Jf^ Uk +>n G S+1 , (£) QL,+,m. (2-55) 
V u R,+,n 

After some manipulation one can obtain a Fisher-Lee expression that generalizes the ID 
expression of Eq. 1.83 

*n,m = ^\j VR '\^ a L R + ' m 4,+,n G S+1 , Q {E) U^. (2.56) 

The notation is a bit messy due to the fact the modes viL/R,+,m/n are n °t orthogonal; see Sec. 
2.2.1. If they were orthogonal, then the expression would simply contain uj^ , n Gs+i ) o(.E)ui, j + )m . 
This indicates that the transmission amplitude from mode W-L,+,m in the left lead to mode u^ + ri 
in the right lead is determined by the Green function matrix that brings you from layer i = 
to layer i = S + 1, i.e. across the scattering region. Proper orthogonalization of the modes 
leads to expressions for , n and u' L + m . The details are a bit messy and can be found in the 
literature. The Bloch velocities f_L/_R 5 + im / n serve to make scattering matrix unitary, as before. 
Since t>L/.R,+,m/n 7^ only for propagating states, Eq. 2.56 expresses explicitly that the trans- 
mission is zero whenever n or m describes an evanescent mode. The layer thicknesses a L / R are 
just normalization factors, since our modes are normalized within a layer. 

After some algebra the Green function expressions we obtained in the one-dimensional case 
can be generalized to three dimensions. The velocity of Eq. 1.77 becomes a velocity matrix 



VL/ fi = ^ImS i/B . (2.57) 

It is a diagonal matrix of dimension N (the total number of modes). The diagonal matrix 
elements are the mode velocities VLm±^ n of Eq. 2.34. Since Vl/b,±,ti = for evanescent states, 
this means that in general the velocity matrices are singular. The transmission amplitudes of 
Eq. 2.56 can be assembled in a transmission matrix t, which using the velocity matrices can be 
expressed as 

t = 2iV-Im£ji G s+ i,o(£) V-ImSi. (2.58) 
Note that this generalizes Eq. 1.84. Finally, the total transmission can be expressed as 



Tr 



ttt 



4Tr [ImEfl G£ +lj0 (£) ImS L Gg )S+1 (£7)] (2.59) 



which generalizes Eq. 1.85. This is known as the Caroli expression or the NEGF expres- 
sion. 25 



25 NEGF = non-equilibrium Green function. In the linear response regime one can use equilirium Green func- 
tions, however, which is what we are doing here. 



